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Abstract 

'a 

' An attempt is made to determine the difference between Bose-Einstein condensation 

of particles and quasiparticles. An equation is obtained for the number of particles in a 
Bose-Einstein condensate as a function of the total number of particles in the system. 
This equation is also written for quasiparticles taking account of their creation by 
external pumping and the presence of equilibrium thermal excitations in the system. 
Analyzing both equations, the chemical potential of the pumped quasiparticles and 
q \ their number in the condensate are found as a function of the pumping intensity. 

A condition under which the Bose-Einstein condensation of low-energy quasiparticle 
excitations starts and occurs at any, including quite high, temperatures is found. 
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1. It is well known that the Bose-Einstein condensation (BEC) is one of the remarkable 
macroscopic quantum phenomena which can occur in one or another collective of particles 
or quasiparticles (QP) with integer spin. There are now a number of examples of the more 
or less direct experimental observation of BEC, among which we call attention to BEC 
in rarifled gases of Bose atoms (Refs. 1 and 2; see also the review Ref. 3). Its critical 
temperature is extremely low (~1CU 6 K) because of the low density of the particles and 
the relatively large mass of even the lightest atoms. In this sense, quasiparticles are more 
interesting. Quite high quasiparticle densities are experimentally fully achievable and the 
effective masses, as a rule, are of the order the electronic masses (to say nothing of the 
existence of massless qausiparticles 4 ). Possibly, tor this reason the question of the BEC of 
QP for the example of large-radius excitons (and biexcitons) was raised about 50 years ago, 5 
and the investigation of the their condensation became very popular especially for indirect 
excitons or in heterostructures. 6 The recent experiments on BEC of excitonic polaritons in 
microcavities also merit equal attention 7,8 . Here the condensation temperature, though it 
increases by an order of magnitude, is still quite low — about several degrees or fractions of 
a degree. 

As concerns the short-radius dipole-active excitons, their condensation has turned out 
to be practically unobservabie because of their short lifetime relative to radiative decay and 
therefore the impossibility of attaining the densities required for BEC. The high-density lim- 
itation substantially weakens for magnons, which also can be regarded as a type of Frenkel 
excitons (or, which is the same thing, short-radius excitons), which possess a long lifetime. 
Apparently, it is this fact that helps to achieve such high magnon densities, so that their 



condensation, as asserted in Ref. 10, is observed in perfect yttrium-iron garnet films with 
temperatures which are anomalously high and compared with room (!) temperature. The 
works mentioned above and' the works of Demokritov's group 11,12 devoted to BEC magnons 
make it necessary to analyze whether or not the conditions at which BEC of QP becomes pos- 
sible at truly high temperatures are indeed satisfied, since nothing like this can be achieved 
for particles and excitons. 

2. First, we note a remarkable feature of BEC, which, being essentially a phase transition, 
can occur in an ideal gas, which makes it possible to calculate many physical charac-teristics 
and observable quantities of this process to a high degree of accuracy. In addition, although 
objects which do not interact with one another at all are unlikely to exist, the concept 
of ideality could be even more characteristic for QP, which interact quite weakly even at 
relatively high densities. 

Nonetheless, it is not ideality (or weak ideality) that lies at the basis of the main difference 
between particles and QP of the Bose type. The fundamental difference between them is 
that if at a given density the particles are characterized by a finite value of the chemical 
potential //, then \i = for equilibrium QP, since their density itself is determined by only 
the temperature. This means that the realization and observation of BEC of QP requires 
converting the system (for example, by special nonthermal action) into an excited (in other 
words, a nonequilibrium) state, which remains in existence for a sufficiently long time (in 
any case, appreciably greater than the equilibration time and the QP lifetime). Then the 
(quasi) equilibrium distribution of QP which is established and which corresponds to BEC 
becomes accessible for experimental study. But, on the other hand, it also follows from this 
that one must take great care when describing the BEC of QP because of the problematic 
nature of using the thermodynamic approach and the nonstrictness of the application of 
various limits. 

We note that there have been more than a few attempts (see, for example, Refs. 13-15) 
at describing the kinetic behavior of the nonequilibrium QP collective (specifically, ferro- 
magnons) under one or another form of pumping — pulsed, noise, and so forth, which can 
then condense. In such strict approaches, additional information must be taken into account: 
on the interaction of QP with one another, their interaction with other objects — phonons 
or defects, the character of the pumping, the temporal evolution, which makes the model 
and the calculations quite complex. All this, being undoubtedly important, nonetheless does 
not appear to be determining for BEC itself. In complete correspondence with Bose and 
Einstein, we understand it as a phenomenon of temperature redistribution of a prescribed 
number of particles over the states and the accumulation of a substantial fraction of them 
(particles) in the lowest of them (states). Even if we confine ourselves to the assumption 
that pumping has created only QP and that fast relaxation made them essentially equilib- 
rium states, we can raise many questions which follow from experiments. 10-12 Neither is our 
objective to examine the coherent manifestations of the QP collective accumulated in this 
manner in the lowest (excited) state, which, certainly, is an interesting and topical problem. 
We shall study the particularities which make it possible for QP to condense at such high 
temperatures, including, in this case, room temperatures; this problem is limited and has 
not yet been thoroughly analyzed. 

3. It is well known that the long- wavelength elementary quasiparticle excitations deter- 
mine the low-energy part of the spectrum of any particular system. On the other hand, their 
interaction with one another is the weakest interaction. As mentioned above, this makes it 
possible to consider, to a first approximation, the subsystem of QP to be an ideal gas. Let 



e(k) = Eo + [e(k) — sq] = eo + £fcj n (k) be the QP energy, and let Eq = e(k ) correspond 
to the lowest excited state. Then, according to the assumption that the number of QP is 
conserved as a result of their external interaction or, in other words, the fact that they have 
a finite chemical potential, we write the number of QP in the kth quantum state at fixed 
temperature as 

n{k ^ ) = e W {[E(k)-^]/k B T}-l- (1) 
Instead of the quantity /i, we introduce, following Ref. 16, a different thermodynamic variable 

H ° [exp(e - /J>)/k B T] - 1' 

whence \i = £ Q — k B T\n[(n + l)/n ]. Using the expression ([2]). the occupation numbers ([T]) 
can be easily rewritten in the form 

1 1 

n(k, fj) -> n(k, n ) = 77—77 s f , T > ^ 7 = n k (T)—- — , (3) 

(1 + I/710) exp[e A . m (k)//c i jrj - 1 1 + [n k (T) + lj/n 

where we have introduced the following notation for the number of particles 

MT) = exp[E kl M/k B T]-V (4) 

which, as one can see from its definition, does not depend on the gap in the QP spectrum. 
The replacement (j2J), eliminating the QP chemical potential from the analysis, facilitates 
the study of the most interesting and important situation where the nonphysical asymptotic 
limit n — > 00 corresponds to the generally accepted condition for BEC fi —> Eq. 1 ^ 

For sufficiently large values of n (3> n k (T)), which, generally speaking, should correspond 
to the BEC regime, the following expansion becomes valid: 

which is a direct indication that the quantity rik(T) introduced above is simply the maxi- 
mum possible (corresponding to the formal condition uq — > 00) occupation number of the 
corresponding state for fixed T. 

It is convenient to use the expansion (J5]) when writing thermodynamically equilibrium 
quantities. Specifically, the total (average) number of particles becomes 

N = J2 9kn(k, no) =n + ShfliK n o) = n + N exc (T, n ), (6) 

where is the degeneracy of the state with wave vector k, and 

n k (T) + l\ , /n k (T) + l\ 2 



N exc (T,n ) = 9kn k {T) I 1 - ( - 
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N exc (T) ~ ^^1, (7) 
n 



corresponds to the total number of QP in the excited states. In the latter equality, the 
following number is singled out: 



N exc (T) = ^ k n k (T), (8) 

which gives the maximum number (achieved under the same condition (see above) n — ► oo) 
of all thermal excitations, and the coefficient 

5N exc (T) = J2 ^ k (T)[n k (T) + 1], (9) 

k^k 

is determined, to first order in l/n , by the fluctuations of the quantity (0). It is evident 
that irrespective of the specific form of the QP dispersion law £fcj„(k) the quantities N exc (T) 
and SN exc (T) are monotonically increasing functions of the temperature, which is easily seen 
from the definitions, since as T increases, all occupation numbers increase. 

Substituting (7) into (jSJ), taking account of © and (jSJ) and assuming N ^ 1 (and, 
correspondingly, n ^ 1), we easily obtain the equation 

N ~ n + N exc {T) - 5Nexc{T) (10) 

n 

for finding the number n of Bose-condensed QP. We find from this equation that the behavior 
of the desired number 

n = n (T) = \{n- N exc (T) + J[N - N exc {T)f + 48N exc (T)} (11) 

for large N and small 5N exc (T)/N 2 changes quite sharply when the temperature crosses a 
certain value T c which is the solution of the equation 

N exc (T c ) = N, (12) 

The Eq. ( II ip shows directly that the critical temperature introduced in this manner depends 
on the total number of particles in the system and the specific form of the function N exc (T). 
The explicit expression for the latter (as also for T c ) is determined by the spectrum, the 
dimensionality, and even the form of the system as well as by the conditions at the boundaries 
of the system. 

If we now introduce the dimensionless density tibec(T) ee n (T)/N of the Bose conden- 
sate, then in the limit 5N exc (T)/N 2 —>■ it becomes a nonanalytic function of the tempera- 
ture: 



n B Ec{T) 



N exc (T) 


e 


N exc (T) 


N 




N 



(13) 



where 9{x) is a step function. Such behavior of the density (fl2|) makes it possible to interpret 
it as the order parameter of the BEC process in a given collective of particles or QP. We shall 
not dwell on this, but we shall focus our attention on the question, raised above, concerning 
the particularities which distinguish QP from particles. 

4. Indeed, thus far the analysis could refer to both of these objects, since the chemical 
potential was excluded from it. However, it is well known and has already been mentioned 



that for QP jj = 0, to that any system at finite T will contain one or another — thermally 
excited — quantity of them (compare Eq. 



^ = - P < 14 > 

In complete analogy with Eq. (J2j) , it is also possible to single out the equilibrium thermal 
occupation of the lowest state: 

4 h \T) ee nt h) = 1 (15) 

exp{so/k B T) - 1 

which is a given (increasing) function of the temperature. 

Using the relation ( TH1) . the numbers ( |T3l) can be easily put into a form analogous to (J3]): 

< h \T) -> 4 h \rMo h) ) = n k (T) - , n * (r)+1 , (16) 
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where, remarkably, the notation ny_(T) is identical to the occupation number (j3J) given above. 
Formally, the expression ffl5|) is exact, but the number of thermal excitations, generally 
speaking, does not have the quantitative advantages over the expression used in the expansion 
(7) (at least, the leading terms) over the occupation numbers n^iT) (or, especially, the 
numbers n^ h \T)) as no over n-^iT) in the BEC regime. However, it can be supposed that 
because the occupation numbers of the excited states decrease exponentially, the error of 
the series expansion in powers of 1/nj) wn l n °t be substantial. 2 - ) Of course, for small 
the difference between the particles and QP becomes less noticeable. However, if one stays 
within the assumption that the initial temperatures are sufficiently high (or energy of the 
elementary excitations is relatively low), then the expansion can be written approximately 
as (compared Eq. (7)) 

TOO - MT) { 1 - + ( ^fi) 2 - ■ • ) , (17) 



and Eq. (fT7|) which follows immediately 



N th ~ 4 h) + N exc (T) ~ (18) 

where, as mentioned, the quantities N exc (T) and 5N exc (T) do not change, even if nonthermal 
(specifically, pumped) QP are present in the system. 

We now introduce their number according to the obvious expression 

n k (T) = nf \T) + n { r mp) (T) (19) 

for each k. Then, according to Eqs. (7) and (jT7|) it is easily shown that for all k ^ k 



nt mp) (T)=n k (T)[n k (T) + l] 



n% h \T) n (T)\ [ n (T) 



1 + n k (T) + 1 



i -i 



\ | n k (T) + 1 
. nf)(T) j 



-i 



(20) 



The latter expression simplifies substantially provided that the proposed (and actually ex- 
perimentally realizable 10-12 ) ratio of the parameters of the system obtains: k B T ^> 6q, fi. 
Indeed, we find from Eq. ( 1201) . taking account of the definition (J2J), that in this case the 
increment to the number of excitations created in all states, with the exception of the lowest, 
by pumping can be represented in the quite simple form 

4™ } (T)^^ k (T)K(T) + l]. (21) 

In other words pumping increases the initial occupation numbers somewhat, but very little 
because of the "temperature suppression factor" that appears \ijk B T <C 1. 

At the same time the density of the pumped QP in a Bose condensate always behaves 
completely differently (see Eqs. (|SJ) and (1151) : 



(„ )(T) = (T) _ n (t h){T) w ksTf^_ (22) 



or can be arbitrarily large because of the possibility /j, — > Eq (except for the exact equality 
H = Eq). It is remarkable that, conversely, in contrast to all n^{T) with k ^ k the occupation 
numbers n < Q Ump \T) increase additionally on account of (compare Eq. (121jl ) the "temperature 
intensification factor" k B T/eo ^> 1. 

5. The role of pumping is obvious — to create a definite number of QP. But how is it 
related with the chemical potential? To answer this question, we take account of the fact 
that (see Eq. (jB])) the total number N of QP in the system can also be divided into two con- 
tributions: the temperature-determined equilibrium part N t h and the pumping-determined 
number N pump : N = N t u + N pump . As is well-known, 18 under electromagnetic pumping 
with intensity I pump this number changes in accordance with the simplest balance equation 
dN/ dt = I pump — (N — N t h) I Trei where r re i is an effective equilibration time after pumping 
is switched off (a more accurate analysis is given in the Appendix). It follows immediately 
from this equation that because of the obvious equality 3 ^ dNth/dt = actually becomes 
dN pump /dt = Ipump — N pump /r re i and that the stationary number (for times t ^> r re i, when 
dN/dt = 0) is the total number of pumped particles, equal to NL^ = I pU mpTrei- Assuming 
that the main number (1221 of the latter accumulates precisely in the lowest state, or that 

T3t 

pump i 



Uq (T) ~ Np U , we easily arrive at 



T - - kBT V 0X\ 

-L pump 'r el / \ > V "J 

E (E - fl) 

which enables us to write the chemical potential of the pumped (and only pumped) QP: 

— P IpumpTrel / Q ,s 

\Kb-L /Eq-\- tpumpTrel) 

This expression shows that for relatively weak pumping, when I pump T re i <C k B T/£ , the chem- 
ical potential of the QP is fi ~ e pU mp T r a / 'k B T ', and for strong pumping when Ip Um p T rei 
k B T/e Q {^> 1) its value goes (see Eq. (jUj) ) to its limit e , always satisfying the physically 
necessary condition fi < e . It is easy to show that the difference 

_ ^ k B T _ k B T ^ k B T 

^0 A* ~ a j s t * (pump) /r r \ ~ J rr ' ^ ' 

ly pump TIq [1 ) J-pump'rel 



and, in consequence, as follows from Eqs. ( 1231) and (1251) . the thermodynamic properties of the 
collective of QP condensed as a result of BEC should depend only on the ratio I pU m P T r ei/kBT. 
It appears that such scaling behavior can be checked experimentally 

Finally, using Eqs. (7) and f|T8|) . and subtracting the latter from the former, gives 

n - n^ h) + — — t^— 5N exc (T) = I pump T re i, (26) 
n n y 

where the unknown is the number Uq(T) of QP in the condensate (or, which is the same thing, 
the number n < Q Urnp \T). Evidently, uq(T) = riQ h \T) in the absence of pumping. In addition, 
we find from Eq. ( l26|) thai, specifically, for weak pumping, when n? wmp) (T)<nf ) (T),the 
number of QP in the lowest level increases as 

no(T) « 4 h) (T) + (l- ^-A_8N exc (T)) I„r rel . (27) 

For strong pumping and when nonequilibrium has been reached n$ ump \T) > nf } (T), the 
number of QP in the lowest level is given by the relation 

n (T) « n£ ump) {T) = I pump r r ei - y^SN^JT). 

k b 1 

Although it is obvious that the latter number also increases as the pump intensity, a "transi- 
tion" occurs from one straight line onto another displaced below the origin of the coordinates, 
which also can be checked experimentally. 

However, if the small corrections to the linear behavior of the number of QP in the 
Bose condensate as a function of the pumping intensity are neglected, which is precisely 
what was assumed in the derivation of the chemical potential, then we arrive at a seemingly 
paradoxical result: the pumped low-frequency QP are always in the BEC regime. However, 
one must remember that the so-called temperatures must satisfy the inequality k B T ^> e . 

This can be stated differently: the computed change of the occupation numbers of the 
quasiparticle states (in the experiments Refs. 10-12 this is a ferromagnetic plate exposed to 
external pumping) shows that the initial — thermally equilibrium — QP "prevent" the accu- 
mulation of new QP in all states except the lowest one (see Eqs. (12T]) and ff22l ). For this 
reason, it seems that pumping-created QP have nothing to do but to collect exclusively in 
it. As follows from our analysis, this will occur at all times (of the order of the action of the 
pumping pulse) irrespective of the temperature as soon as, once again, k B T ^> £ - Thus, 
the satisfaction of this inequality makes the temperature a factor which not only does not 
impede (neglecting relaxation processes, which decrease r re /) but rather even promotes BEC 
to some extent. Of course, as mentioned above, the latter assertion relies on the assumption 
that the temperature remains constant. In addition, temperature stability in a real experi- 
ment requires special attention. However, it is worth repeating that when high values of the 
temperature are reached, its relative increase as a result of irradiation processes should not 
be too large. 

6. The results obtained above show that certain types of QP, specifically, QP with low ac- 
tivation energy (or at temperatures not exceeding them) under external pumping can (more 
accurately, must) accumulate only in the lowest possible state, which corresponds completely 
with the BEC phenomenon. In addition, such accumulation is accompanied by relaxation 



from the strongly nonequilibrium state into (thermal) a (quasi) equilibrium state which is 
found to be indistinguishable from a Bose-condensed state. It appears that experimental 
verification (including by measuring the Mandel'shtam-Brillouin scattering 10-12,19 ) of the 
computational results obtained in the present work could show that there is a large differ- 
ence between BEC particles (for example, Bose atoms of light metals) and low-frequency 
QP. This difference could be manifested, first and foremost, even in its (BEC) most im- 
portant characteristic — the critical temperature. Other features (specifically, coherent and 
fluctuation properties) of BEC for QP of different nature, with different dispersion laws, and 
in systems with different dimensionality will be analyzed elsewhere. 
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APPENDIX 



To simplify the exposition the balance equation given above for determining the number 
of pumped QP was written in a form which gives only the proportionality between A 7 ^ and 
the absorbed power I pU mp T rei] there are no terms which give rise to a pump threshold near 

which N pump ~ ^Jlpump- 12,20 

In this connection we recall that under the conditions of parametric electromagnetic 
excitation the exponential growth of the number of magnons in the system occurs when the 
creation rate of these QP is higher than their annihilation rate, whose determining channel 
are pair collisions. 21 Then, to take these circumstances into account we shall write the balance 
equation in the form (see Ref. 19) 

**L-T C\ T l "O N ~ N * ( N ~ N ^) {2) 

, J-pumpK 1 * T L) /jn 2 



T rel rel 



or, in accordance with the notation adopted, 



2 



dNp Urnp _ ( l\ ^ N pW)lj) 

^ J-pump ~r \ -ipump ^ I ly pump ' V^ 11 / 

\ ^rel ' T re i 

where T pum p = Ipum P (N t h + 1), and t^J and t^J are the lifetimes of the magnons with respect 
to one- and two-particle processes. 

It is easily shown from Eq. (Al) that the threshold behavior of the number of pumped 
QP is completely determined by the sign of the second term on the right-hand side of this 
equation, and the last term limits the growth of their number. It is precisely the equality 
Ipump = V T rli that sets the magnitude of the desired threshold, as happens in the S theory 
of parametric resonance. 21 

As for the stationary number of pumped magnons, according to the same equation (Al) 
^pump ~ Ipumpi~rl} which obtained above; near threshold, when I pU mpTrei ~ 1> ^ ne be- 
havior changes and Nf ump ~ \jlpumprfl ] finally, if I pump > l/r r ( ^, the number Nf ump ps 
( 2 .) lo j- . Frr J?) Jo\2 _l T J?) 



IpumpT re i/2 + J (IpumpTrJ /2) 2 + I pum pT re \ , i.e. it is described by a dependence which is, once 



again, close to linear. This is the qualitative picture observed experimentally, where the 
function (I pump ) ~ n^ mp (I pump ) changes from linear to square-root back to linear 

behavior. 20 If necessary, the corresponding dependences can be easily introduced into the. 
relations fl23|) - fl26|) . whose meaning remains intact. 

a )E-mail: abugrij@bitp.kiev.ua 
6 ^E-mail: vloktek@bitp.kiev.ua 

1 - ) We note that in the theoretical study in Ref. 17 of the fluctuations of a Bose-condensate 
of ultracold atoms the chemical potential was also eliminated from the expression for the 
number of particles but in a different manner so that the noted divergence of the quantity 
remained. 

2 ^Here we implicitly assume that ^> 1, which can hold only under the obvious con- 
dition ksT ^> 6q and holds with room to spare, for example, in experiments with BEC of 
ferromagnets. 10-12 At the same time ~ exp(— e/fceT) <C 1 almost always holds for ex- 
citons or QP with a large gap, whence it follows that the Boltzmann distribution describes 
their thermal number. Here n£ h (T) « n^{T) exp(— So/hsT) <C n^{T). 

3 **0f course, neglecting a possible change of temperature, whose relative magnitude near 
10 2 K cannot be large. 
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